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Abstract 

The movement of the electrons under the simultaneous influence of 
a scalar periodic potential and of a uniform transversal magnetic field is 
described by the well-known second order discrete Harper equation. This 
equation originates from a two-dimensional energy dispersion law under 
the minimal substitution. Here one deals with the Harper model under 
the additional influence of the discretized spin orbit interaction. Convert- 
ing the spin-orbit interaction in terms of discrete derivatives opens the 
way for analytical and numerical studies. One finds coupled equations 
for the spin dependent wave functions, which leads to an appreciable 
alteration of the nested energy subbands characterizing the self-similar 
structure of the usual Harper spectrum. To this aim the transfer matrix 
method has been applied to selected spin-up and spin-down wavefunc- 
tions. Accordingly very manifestations of spinfiltering and of spin corre- 
lations are accounted for. Our energy-bands calculations show that the 
splitting effect implemented by such wavefunctions is appreciable. 
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1 INTRODUCTION 

Two-dimensional systems of electrons submitted to the simultaneous in- 
fluences of the periodic potential and of a homogeneous external magnetic field 
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have received great interest [TJ EJ in the last four decades. Such studies led to 
the celebrated Harper equation PQ: 



V-Wi + tPn-i + 2Acos(27r— n + k 2 a)t(; n = E s tp, 



n 



(1) 



which has been extensively investigated. This equation originates from the 
influence of the minimal subtitution on the energy dispersion law. Here E s = 
2E/Eq denotes the dimensionless energy, which comes from a single band tight 
binding description with nearest neighbor hopping. 

Further, A denotes the anisotropy parameter, while (3 = (j)/4>Q stands for 
a commensurability parameter expressing the number of flux quanta per unit 
cell. The magnetic flux and the flux quantum are denoted by 4> = Ba 2 and 
4>q = h/e respectively, while k 2 a stands for the Brillouin phase. 

Equation (pQ), is a second order discrete equation which has been studied 
numerically in terms of the transfer matrix method by Hofstadter [3]. It has 
been found that the spectrum is characterized by a nested band structure i.e. 
by the so called Hofstadter-butterfly. Besides experimental relevance [I], the 
Harper-equation provides ideas for useful applications concerning electrons on 
lattices under the influence of a transversal magnetic fields (0,0, B) [5]. 

On the other hand, it has been demonstrated that the influence of the 
spin-orbit interaction (SOI) E] produce sensible effects [8]. Such results 
open the way of controlling the electron's spin under the influence of an ex- 
ternal electric field or of a gate voltage, too [9]. 



In order to investigate the influence of the Rashba-SOI (RSOI) on the 
Harper model, we shall supplement the Harper Hamiltonian with the dis- 
cretized version of the Rashba energy term: 



where A is the Rashba coupling parameter, a x and a y are the Pauli matrices, 
while p x and p y stand for momentum operators, as usual. 

This later discretization has its own interest as it provides the possibility 
to account for the influence of an underlying electric field in a rather consistent 
manner. This differs, of course, from the magneto-electric Harper equation 
discussed before [10j[Tll[12]. I n this paper one deals with the numerical energy- 
bands realization of the superposition of (1) and (2), now by accounting for 
selected wave functions. Such solutions concern spin- up and spin-down wave 
functions, respectively, but inter-related ones will also be considered. 

The paper is organized as follows. The discretized RSOI model as well as 
the derivation of the coupled equations are presented in section II. Numerical 
investigations of the coupled equations are done in section III by resorting 
to some selected wave functions. First one accounts solely for the influence of 
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the spin- up and spin- down wave functions, respectively, which is reminescent 
to spin-filtered systems. Next, one deals with inter-related spin-up and spin- 
down wave functions. The conclusions are presented in section IV. 



2 DERIVATION OF COUPLED EQUATIONS 

In order to implement the discretized version of RSOI into the Harper 
model, we have to resort to the modified interaction Hamiltonian 

~ _ . / b\ - id 2 + eBx/h\ . . 

M °° ~ CA [- dl _ i d2 + eBx/h J [6) 

where e > 0, d\ = d/dx and d 2 = d/dy. For this purpose one proceeds by 
applying the minimal substitution to by using the Landau-gauge for which 
the vector potential is given by (0, Bx, 0). 
Next, one resorts to the factorization 

^x,y)=eMihy)(^ { $) (4) 

where ip^(x) and ipi(x) stands for spin-up and spin-down wave functions, re- 
spectively. This leads to the iiS'O-Hamiltonian 

where 

L>i ±) = ±dt + k 2 + ^Bx . (6) 

On the other hand, the Harper-equation is generated by the energy dis- 
persion law 

Edisp(ki,k2) = -Eb(cos k\a + A cos k 2 a) (7) 

under the influence of the Peierls substitution. This results in the substitution 
rules k\ — ► —id\ and k 2 — » — id 2 + eBx/Ti. These rules are responsible for the 
onset of the Harper-Hamiltonian such as given by 

H H ar = : y[2cosh(a9i) +2Acos(^^ + k 2 )} . (8) 

One would than obtain the second-order discrete equation 

2 ~ ae 

—H Ha ri>-[A{x) = ip^ l (x + a) + il;^ l (x-a) + 2Acos(—Bx + k 2 a)ip^ l (x) (9) 
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which reproduces equation ([T]) in terms of the space-discretization x = na, 
where n is an arbitrary integer. What then remains is to establish the dis- 
cretized version of equation ([5]), by proceeding via 

d 1 f(x) = 5 1 f{x) + 0{a 2 ) (10) 

where 5\ stands for discrete derivative. It is understood that 

<5i/(x) = - sinh(a9 1 )/(x) = ±-(f{x + a) - f(x - a)) . (11) 
a la 

Accordingly, d\ and <5ibecome identical to first a-order. This shows that 
the discretization one looks for, is provided by the modified substitution rule 

£>j ±) ->±$ 1 + Jfcjj + 27rn-— (12) 

a 0o 

which stands for the discrete counterpart of equation The eigenvalue 

equation characterizing the total Hamiltonian is then given by 

Statist, j. 0*0 = {Huar + H s0 )ip^ |(x) = Eip 1t (13) 

where E denotes the total energy and where H so stands for the discretized 
-RS'O-Hamiltonian. This results in the coupled equations 

E s ip^(n) = </3|(n + 1) + (/?|(n — 1) + 2A cos(27r/3n + /c2a)(/?|(n) + 

r r 

+— [ip^n + 1) - ip^n - 1)] + ip l (n)-(k 2 a + 2tt f3n) (14) 

with 

E s ip^(n) = </9|(n + 1) + (/?|(n — 1) + 2A cos(27r/3n + A;2«) ( P|(n) + 

r r 

+ — [-9J T (n + 1) + if^{n - 1)] + ip 1 (n)-(k 2 a + 2vr/3n) (15) 

where T = —2A/Eq contains the Rashba coupling parameter. In order to 
highlight the meaning of the coupled equations, we shall resort, this time, to 
some quickly tractable simplifications. 

3 NUMERICAL INVESTIGATIONS 

At this point we have to realize that Eq can be established in terms of the 
effective mass m* of the electron as Eq = —Ti 2 /m*a 2 . Using m* = 0.067m e and 
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a = lOnm [13] then gives Eg ~ 11.37meV. In addition one has A ~ 5-10 -11 eVm 
[14:\ [T5] . A first simplification is to assume that spin-down wave functions are 
zero (Pi(n) = 0, which can be viewed as relying on spin filter devices [16j or 
spin separation processes [17] . 

This Ansatz leads to a s = 1 spin-selection and # 2 = k<ia, as 

<p^(n + 1) + </2 T (n - 1) = -2Acos(2vr/5n + 2 ) + E s ip^{n) (16) 

with 

<^ T (n + 1) - p T (n - 1) = 2(2vr/?n + 2 )</? T (ra) (17) 
which proceed irrespective of T. Under such conditions we obtain 

¥> T (n + l) =y> T (n)[^-AcoB(27r/3n + fl 2 ) + (27r/9n + e 2 )] (18) 



with 



so that 



¥> T (n-l) =(^ T (n)[^-Acos(2^/3n + 2 )-(27r/3n + 2 )] (19) 
V? T (n - 1) _ £i s - 2A cos(2vr ( 5n + 2 ) - 2(2vr/?n + 2 ) 



and 



<£ T (n) 
y? T (re) 



This amounts to consider 



where 



^- „ 1 „ ^ 



(20) 



(21) 



V? T (n + 1) £ s - 2A cos(2vr/3n + 6> 2 ) + 2(2vr/3n + 2 ) 
Applying the standard transfer matrix method, one obtains 

r m(oh and ^T(2)) =TT (n(oh (22) 
which can be generalized via 

(^^^■••^O- (23) 



/ i^ T (n) \ rp (<Pf(n-l)\ (l)A \ 

Wn+lV ~ ±n \ c^(W> > ^> 



T = 2 -Acos(27r / 3n+e 2 )+(27r / 3n+e 2 ) 

' '" q g s -2A cos(27r/3n+27r/3+6> 2 )--2(27r/3n+27r/3+e2) i 

(25) 



5 



supposing that 

\Tr{T n )\ = |2Acos(27r/?n + 6> 2 )| < 2A . (26) 

Note that A = 1 case, which will be assumed hereafter, represents the 
critical point of a metal-insulator transition. We are by now in position to 
perform numerical studies in accord with (I24|) -(j25"j). One realizes that this 
time, the usual symmetries of the Harper-spectrum like such as the exact 
energy reflection symmetry E — > —E and the /3-symmetry, i.e. the symmetry 
under (5 = 1/2 — > (3 = 1 — 1/u (u = 1, 2, 3, 4, ...), cease to be fulfilled. 

The understanding is that by virtue of the s = 1-selection one deals with 
a different experimental situation. The flux dependence of the energy bands 
when s = 1 (i.e. 7 = 0) is displayed in Fig.l, as illustrated by the bands 
located in the upper half plane. However, we can also select tp+ = (instead 
of (p^ = 0), in which case one deals with s = — 1 (instead of s = 1). 



Fig.l The flux dependence of energy bands for s = 1 (upper 
half plane) and s = — 1 (lower half plane). 



This time the energy-bands get located inside the lower half plane. 

One remarks that bands become increasingly dispersed as the magnetix 
flux increeases. This means that (3- symmetry is completely lost. When the 
flux approaches unity, the bands become equidistant, in accord with the be- 
havior of Landau levels. Such patterns can be gathered together as 

2^r = sF ± ((3,k) (27) 
Eo 

in which s = 1 or s = - 1, respectively. In addition there is F±(k) > and 

F+(f3,k)nF„(P,k) (28) 
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which is responsible for energy reflection symmetry. The next approximation 
is to consider inter-related spin- up and spin-down wave functions say V 9 |( n ) = 
j(p^(n), where 7 is a real parameter. 



Fig.2 The flux dependence of energy bands for subunitary values of the 
7-parameter, like 7 = ±1/4, ±1/2, ±17/20 and 7 ± 1. 

One would then obtain 

[E a - 2A cos(2vr/3n + k 2 o) - 2jT*(k 2 a + 2vr/3n)]^ T (n) = 

= (1 + 7 r> T (n + 1) + (1 - lT*)^{n - 1) (29) 

with 

[E s -f - 27A cos(27r/3n + k 2 a) - 2T*(k 2 a + 2vr/3n)]^ T (n) = 

7 



= ( 7 - r> T (n + 1) + ( 7 + r> T (n - 1) 

by virtue of (JT3J) and (fT5|) . respectively. 



Fig.3 The flux dependence of energy bands for 
values of supraunitary modulus of the 7-parameter, 
like 7 = ±1.05, ±1.2, ±2 and 7 = ±99. 



These later equations lead in turn to 

^{n-l) _ (7-r*)^(n) - (l± 7 r*)£ 2 (n) 



p T (n) (7 - r*)(i - 7 r*) - (1 + 7 r*)(7 ± r*) 

and 

(7 + r*)(i + 7 r*) - (1 - 7r*)(7 - r*) 



¥» T (n+l) (7 + T*)£?i(n) - (l-7r*)E 2 (n) 



where F* = F /2a so that F* ~ 7.043 • 10 20 is vanishingly small. In addition 



Exin) =E a - 2Acos(2vr/?n + k 2 a) - 2jF* (k 2 a + 2?r/3n) 



(33) 



and 



E 2 (n) = E s ~/ - 27A cos(2vr/3n + k 2 a) - 2F*(k 2 a + 2ir(3n) 



(34) 



denotes the energies from (I3ip and (|32|) . Using again the method of the transfer 
matrix, one finds energy-band realization which are rahter sensitive to 7. In 
addition we have to remark that by now E s = 2E/Eq is an even function of 
7. When I7I < 1, the energy bands exhibit dispersions decreasing with |7|, 
as shown in Fig.2 for 7 = ±1/4, ±1/2, ±17/20 and 7 = ±1. We have to 
remark that in this latter case energy dispersion effect get inhibited, now by 
preserving both energy reflection and /3-symmetries. We have also to remark 
that the plots displayed in Fig.2 for 7 = ±1 reflect, within a reasonable 
approximation, the onset of a two level configuration, say 



proceeding irrespective of the magnetic flux. Such configurations are of a 
special interest for applications in the field of optoelectronics |18} I19j. Next 
we found that supraunitary values of the 7 parameter lead to an asymmetric 
location in the flux dependence of energy bands, as shown in Fig.2, for 7 = 
±1.05, ±1.2, ±2 and 7 ± 99. One sees that in the later two cases the energies 
exhibit positive values, only. 



In this paper we have investigated both theoretically and numerically the 
competition effects between the Harper model and the discretized RSOL 

Coupled equations for the spin dependent wave functions have been es- 
tablished in some particular cases. Unlike the Harper spectrum, which, is char- 
acterized by exact symmetries mentioned before, the presence of the RSOI 
leads to dispersion effects in the flux dependence of the energy bands, except- 
ing the 7 = ±1 plots in Fig.2. 

Using the transfer matrix technique, we found that the RSOI produces 
dispersive effects in the energy dependence on the magnetic flux. Such effects 
get visualized by a sensible alteration of the rather symmetrical nested energy 
sub-bands characterizing the self-similar structure of the usual Harper spec- 
trum. We then have to realize that, excepting of course the special 7 = ±1 
case mentioned above, both energy-reflection and /3-symmetries are lost. Our 




±2.36 



(35) 



4 CONCLUSION 
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first but restrictive approximation is to account only for the influence of spin- 
down and spin-up wave functions, respectively. The modified energy bands 
established in this manner are displayed in Fig.l. 

Our next approximation is to account for inter-related spin- up and spin- 
down wave functions. Now we found that energy bands exhibit again visible 
modifications, as shown in Fig.2 and Fig.3. 

In other words we found that supplementing the Harper-Hamiltonian 
with the discretized RSOI leads to significant modifications of energy bands, 
which serves to a better understanding of the electronic structure. This opens 
the way to further applications concerning the influence of the spin-orbit effect, 
with a special emphasis of two-level realizations mentioned above which are of 
interest in the field of optoelectronics. More general wave functions can also 
be considered, but such goals go beyond the immediate scope of this paper. 
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